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ABSTRACT 



Fireballs of gamma-ray bursts are partially made of free neutrons. Their presence 
crucially changes the mechanism of the fireball deceleration by an external medium. At 
late stages of the explosion, neutrons fully decouple from ions and coast with a constant 
Lorentz factor T„. As the ion fireball decelerates, the neutrons form a leading front. This 
front gradually decays, leaving behind a trail of decay products mixed with the ambient 
medium. The kinetic energy of the decay products far exceeds the medium rest-mass 
energy, and the trail has a Lorentz factor 7 3> 1 at radii R < Rtmi\ ~ WRg ~ 10 17 cm, 
where Rg ~ 10 16 (T n /300) cm is the mean radius of neutron decay. The ion fireball 
sweeps up the trail material and drives a shock wave in it. Thus, observed afterglow 
emission is produced in the neutron trail. It can naturally re-brighten or display a 
spectral transition at R ~ -Rtraii where the impact of neutrons turns off. Absence of any 
neutron signatures would point to an extremely low baryon loading of the fireballs and 
a strong dominance of a Poynting flux. 

Subject headings: Cosmology: miscellaneous — gamma-rays: bursts — radiation mech- 
anisms: nonthermal — shock waves 



1. Introduction 

The presence of neutrons in gamma-ray bursts (GRBs) and their possible role for the observed 
emission was realized in the past few years (Derishev, Kocharovsky, & Kocharovsky 1999a, b; Bahcall 
& Meszaros 2000; Meszaros & Rees 2000; Fuller, Pruet, h Abazajian 2000; Pruet & Dalai 2002). In 
an accompanying paper, we study in detail the nuclear composition of GRB fireballs and show that 
the presence of neutrons among the ejected baryons is practically inevitable (Beloborodov 2003, 
Paper 1). One implication is an observable multi-GeV neutrino emission from inelastic neutron- ion 
collisions in the fireball. Here, we focus on a different aspect. We show that the neutrons have 
a dramatic impact on the explosion dynamics at radii as large as 10 17 cm and propose a novel 
mechanism for GRB afterglow emission. 



1 Canadian Institute for Theoretical Astrophysics, 60 St. George Street, Toronto, ON M5S 3H8, Canada 

2 Physics Department, Columbia University, 538 West 120th Street New York, NY 10027 

3 Astro-Space Center of Lebedev Physical Institute, Profsojuznaja 84/32, Moscow 117810, Russia 



- 2 - 



Let us remind what happens in a relativistic explosion without neutrons (see Meszaros 2002 
for a review). The ejected fireball with mass M e j and Lorentz factor T c j sweeps up an ambient 
medium and gradually dissipates its kinetic energy. The dissipation rate peaks at a characteristic 
"deceleration" radius i?dcc where half of the initial energy is dissipated. This radius corresponds to 
swept-up mass m^ cc = M c j/T e j. Further dynamics is described by the self-similar blast wave model 
of Blandford & McKee (1976). How does this picture change in the presence of neutrons? 



2. Neutron front and its trail 

Neutrons develop a Lorentz factor T n = 10 2 — 10 3 at the very beginning of the explosion when 
the fireball is accelerated by radiation pressure (Derishev et al. 1999b). They are collisionally 
coupled to the ions in the early dense fireball, and decouple close to the end of the acceleration 
stage. Then the neutrons coast and gradually decay with a mean lifetime Tp s=s 900 s and a mean 
decay radius Rp = crpT n , 

RfS = 0-8 x 10 16 (J^j cm. (1) 

At radii under consideration, R > 10 15 cm, the ejected fireball is a shell of thickness A <C R. 
In contrast to neutrons, the ion component of the fireball is aware of the external medium and its 
Lorentz factor T decreases. As T decreases below r n , the ions fall behind and separate from the 
neutrons. Thus the fireball splits into two relativistic shells which we name N (neutrons) and I 
(ions); the N-shell will also be called "neutron front". The I-shell lags behind by a distance I, 

2f^~2f2' (2) 

where and (3 n are the shell velocities in units of c. For simplicity, we hereafter focus on the 
stage of complete separation I > A (it sets in right after the beginning of the I-shell deceleration if 
A ~ 10 11 cm, and it covers the whole explosion if A — > 0). 

The mass of the N-shell gradually decreases because of the /3-decay, 

M„(i?) = M n °exp(-A). (3) 

However, even at R > Rp, the N-shell energy E n = T n M n c 2 is huge compared to the ambient rest 
mass mc 2 . For example, at R = R^ecj E n /md ec c 2 = X n T n T e j exp(-i? dcc /i? /3 ) where X n = M°/M ej . 

The neutron decay products p and e~ share immediately their huge momentum with ambient 
particles due to two-stream instability (the instability timescale is set by the ion plasma frequency 
uj p and it is the shortest timescale in the problem). Thus, the N-shell leaves behind a mixed trail 
with a relativistic bulk velocity (3 < (3 n which we calculate now. 

Let dm = (dm/dR)dR be ambient mass overtaken by the N-shell as it passes dR and dM n = 
(M n /Rp)dR be mass of decayed neutrons. The dm and dM n share momentum and form a trail 
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element with proper mass dm* = dM n + dm + dmh ea t which includes heat dissipated in the inelastic 
dm-dM n collision. The laws of energy and momentum conservation read 

T n dM n + dm = 7 dm*, (3 n T n dM n = /3 7 dm*, (4) 

where 7 = (1 — /3 2 ) -1 / 2 is the trail Lorentz factor. Let us denote 

dM n M n (dmy 1 
aR) = ^=R- p [dR) ■ (5) 

From equations (4) we find 

It gives 7 > 1 until essentially all neutrons have decayed. For illustration, let us specialize to a 
power-law density profile of the ambient medium. Then 

R\ k M„ ( R N L ~ A 



and £ evolves from ( = M°e _1 jkmp S> 1 at R = Rp to ( <C 1 at R 3> Rp as M n decays 
exponentially. There exists a characteristic radius -Rtraii where the trail becomes nonrelativistic 
((3 = 0.5). It is defined by condition ( = T^ 1 (see eq. 6), which requires about 10 e-folds of the 
decay (for a typical mp ~ m^ cc ~ 10~ 5 M°). Thus, 



i? trai i « lOi^ = 0.8 x 10 17 



cm. (8) 



-Rtraii depends very weakly (logarithmically) on the ambient density and X n . 

Now that we know (3, the trail density n can be easily calculated. The neutron front has a small 
thickness A and the ambient particles cross it almost instantaneously (on timescale 7 2 A/c <C R/c). 
Measured in the N-shell frame, the flux of ambient particles is (3 n T n rio = (f5 n — /?)r„7n am b where no 
and ra am b are proper densities of the ambient particles ahead and behind the N-shell, respectively. 
This gives the medium compression as it is accelerated in the neutron front: n am b/no = [Pn/liPn — 

The total density of the trail includes the neutron decay products, so n = (1 + C)^amb- Using 
equation (6), we get 

The energy dissipated in the neutron front is given by 

dE n . { . = 7 (dm* - dm - dM n ) c 2 = (T n - 7) dM„c 2 - (7 - l)dm c 2 . (10) 
After simple algebra (using eqs. [4]), we find 



dE, 



n.f. 



FnPn ~ 7/3 

dR \T nl {(3 n -(3) 



dm 2 
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It simplifies to dE Q f = 2^ 2 dmc 2 for the interesting regime 1< 7 < T n (I^ 1 < ( < T n ). 

The trail is formed very hot. This can be seen when comparing inertial mass dm* (that includes 
heat) with rest mass dMj, = dm + dM n = (l + ()dm, which gives dimensionless relativistic enthalpy 

dm, (C 2 + 2r„C+ 1) 1/2 , , 

" = dM = i+c ■ (12) 

For T^ 1 < ( < T n we find /i>l, i.e. internal energy of the trail far exceeds its rest-mass energy. 
The trail parameters are summarized in Table 1. 

It is instructive to view the dissipation process in the rest frame of the trail. Here, the 
elements dm and dM n have initial Lorentz factors 7 and 7 = IVy(1 — /3(3 n ), share their opposite 
momenta, and come at rest. This is achieved via a plasma instability that isotropizes the particle 
momentum distribution. It may result in two isotropic ion populations in the trail: ambient ions 
with mean Lorentz factor 7 and /3-decay protons with mean Lorentz factor 7. Both populations 
are relativistically hot, 7 » 1 and 7 > 1, as long as T" 1 < ( < F n . They have equal energies 
jdM n = 7dm (which correspond to equal relativistic bulk momenta of dM n and dm before they 
come at rest in the trail frame). The ratio of their temperatures is the reciprocal of their density 
ratio. 

A detailed model of dissipation in the neutron front is an interesting plasma physics prob- 
lem, which we defer to a future work. We emphasize here that it is different from dissipation in 
collisionless shocks. The thickness of a shock, S ~ c/u p , is set by the timescale of two-stream 
instability, and it is a discontinuity in the hydrodynamical sense. By contrast, the N-shell has 
thickness A ~ 10 11 - 10 12 cm which is orders of magnitude larger than 5. The neutrons decay and 
cause dissipation everywhere in the N-shell. The ion medium velocity grows smoothly from at the 
leading edge of the neutron front to (3 at its trailing edge. This can be called "volume dissipation" . 

3. Shock wave 

The ion fireball follows the neutron front and collects the trail. As a result, (1) the ion Lorentz 
factor T decreases and (2) a shock wave propagates in the trail material. The shock has a Lorentz 



Table 1: Trail parameters 





C>r„ 


i< C < r„ 


r-' < c < 1 


C<rv 


7 


r„ 


(r„C/2) 1/2 


(r„C/2) 1/2 


1 


n/n 


C 2 


2 7 c 


27 


1 


A* 


1 


r„/ 7 ~ 27 


27 


1 
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factor T ^ T s h < r n and cannot catch up with the neutron front (unless no[R] falls off steeper than 
R~ 3 ). The shock lags behind the neutron front by distance Z s h ~ R{(3 n — fish) an d its time lag is 

tsh ~ c(p sh -0)~ (r s \ - 7 2 ) c ■ (13) 

It is much shorter than R/c as long as T s h 2> 7. Hereafter we use the approximation < R/c 
and assume that the trail is quickly picked up by the shock, before it could expand and change 
its density or velocity. We allow, however, for rapid radiative losses of the trail: it may cool on a 
timescale i coo i <C R/c and possibly t coo \ < i s h- Below we consider two extreme regimes: t coo \ 3> t s h 
(adiabatic) and i coo i <C i s h (radiative with fi = 1). 

In the radiative regime, the energy dissipated by the shock is 

^-rCr--D<i + 0^, <"> 

where r re i = T7(l — f3ij3) is the trail Lorentz factor with respect to the I-shell. For T" 1 <C ( <C T n 
one can use approximate expressions ( = 2^/ 2 /T n and r rc i = (l/2)(7/T + r/7) to get dE^ = 
(I/27 + 7/r n )(r — 7) 2 dmc 2 . Note that the dissipation is smaller than it would be in the absence 
of the neutron front, dE s h = T(T — l)dmc 2 . Hence, in the radiative regime, the neutrons delay the 
deceleration of the ion fireball. 

In the adiabatic regime, the preshock material has a high enthalpy /j, given by equation (12). 
The postshock heat, measured in the lab frame, is d-Eheat = r(r re i/x — 1)(1 + Qdm. It includes heat 
deposited into the preshock medium by the neutron front (eq. 10), which should be subtracted to 
get the energy dissipated in the shock itself. Hence, 

= [r(r rcl ^ - 1) - 7 (/i - 1)] (1 + 0^ c 2 . (15) 

For T^ 1 < C < r n we get dE n L = 2<y 2 dmc 2 , dE sh = (T 2 - 7 2 )dmc 2 , and the total dissipated 
energy d-Ehcat = (F 2 + 7 2 )dmc 2 . As long as T > 7, the bulk of energy is dissipated in the shock 
rather than in the neutron front. Moreover, the dissipated energy is the same as in the absence of 
a neutron front, and hence the fireball deceleration is the same. 

In both radiative and adiabatic regimes, the shock dissipation is suppressed when the blast 
wave decelerates to T ~ 7. At T = 7 the shock would disappear. Thus, T is bound from below by 
7- 



4. Numerical examples 

Let us consider the simple radiative regime. The I-shell and the shocked part of the trail can 
be treated as a single shell with a growing mass M(R) = M c j + m — M n (a known function of 
radius) and a decreasing V. Mass gain dM causes deceleration dr that is found from the energy 
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and momentum conservation: d(TM) = -fdM - dE iad /c 2 and d(ftrM) = (3-fdM - (dE rad /c 2 )(3i 
where dE Ta ^ = dE s h is the radiated energy. Excluding dE ra a, we get the dynamic equation, 

M dm = ~ v2im ~ m+0 ^n- (16) 

Using equations (6) and (7) we can solve numerically equation (16) for T(R) with an initial condition 
r e j. The results depend on the type of the ambient medium. Figure 1 shows two examples: a 
standard ISM and a wind from a Wolf-Rayet progenitor. 

The radiative explosion has two separate emission regions: behind the neutron front and 
behind the shock front. Their luminosities equal the corresponding dissipation rates (cdE Q j /dR) 
and (cdE s h/dR) (Fig. 1). The neutron front dissipation peaks at Rg. The shock dissipation can 
have two peaks (if mp 3> M e j/r e j as in the wind example in Fig. 1). The first peak marks the 
beginning of the I-shell deceleration, and it is followed by a minimum when T approaches 7. At 
R > 2Rp, 7 falls down steeply and the shock becomes strong again (high r re i), leading to the second 
peak at R ss i?traii- 



5. Discussion 

Dynamics of neutron-fed explosions is a clean physical problem: existence of neutrons in GRBs 
and their mean lifetime of 15 min is all we needed to construct the model. We focused here on the 
simple case where all neutrons have equal Lorentz factor r n , and data may require a multi-shell 
picture with variable T n . This extension could affect our model in case of high variations Ar n > r„ 
(though at large R, where neutrons die out exponentially, it is sufficient to consider a highest T™ ax as 
slower neutrons have shorter lifetimes and their population is smaller by exp[— r™ ax /r n ]). Another 
simple extension is beamed ejecta. Then at late stages the ion fireball spreads laterally while the 
neutron beaming remains constant. 

In contrast to dynamics, emission is really complicated. In general, radiation from collisionless 
blast waves is not derived from first principles, and neutrons do not make the problem simpler. The 
afterglow emission is believed to be synchrotron, and it depends on poorly understood generation 
of magnetic field and electron acceleration. The standard model without neutrons relies on the 
field generation in the shock by the two-stream instability (Sagdeev 1966, Medvedev & Loeb 1999, 
Gruzinov 2001). The shock front is, however, extremely thin (S ~ c/uj p <C R/T 2 ), and the postshock 
field decays quickly. The model needs a significant remnant field in an extended layer behind the 
shock, which is uncertain. This problem is alleviated in neutron-fed explosions. Here, the leading 
neutron front is an additional dissipation region maintained in a turbulent state with generated 
magnetic fields. The N-shell has thickness A comparable to R/T^. It can produce a significant 
synchrotron radiation by itself 4 and leave behind remnant fields and hot plasma to be used by the 



4 The N-shell emission decays exponentially on observed timescale RpjcV^ — TjjjY n ~ 1 — 10 s and it can be 
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Fig. 1. — Radiative model with X n = 0.5 and T n = T e j = 300. The ambient medium is parametrized 
by rrip and k (eq. 7). Left panels show the case of k = 1 and mp = 10M e j/r e j (typical for a wind 
from a Wolf-Rayet progenitor). Right panels show the case of a normal interstellar medium with 
k = 3 and mp = M e j/r e j. Top: Lorentz factors of the neutron trail (7) and ion fireball (T). Dotted 
curve shows the fireball deceleration T(R) that would take place without neutrons. Bottom: Radial 
distributions of dissipation rates in the neutron front and the shock. E e j = r c jM c jC 2 is an initial 
total energy of the ejecta. R is measured in units of the mean decay radius (eq. 1). 
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ensuing shock wave for afterglow emission. We emphasize differences from a customary external 
shock: the neutron-trail shock propagates in a relativistically moving, dense, hot, and possibly 
magnetized medium. The fate of magnetic fields in such shocks and the resulting emission is an 
interesting issue for a future study. 

The neutron impact ceases at i?traii ~ 10 17 cm, which can leave an imprint on the observed 
afterglow. For example, the shock dissipation can have a bump (Fig. 1), and a spectral transition 
is also possible. The arrival time of radiation emitted at -Rtraii is ~ -Rtraii/2r 2 c (counted from the 
arrival of first 7-rays). It may be as long as 30 days or as short as a few seconds, depending on the 
fireball Lorentz factor r(i? trai i). Remarkably, -Rtraii is almost independent of the ambient medium, 
and its observational signature would give information on r(i? tra ji). Recent early observation of 
a GRB afterglow (GRB 021004) discovered an interesting re-brightening at 10 3 s which would 
correspond to r(i? tra ii) ~ 30(r„/300) 1 / 2 . Also, we do not exclude a possible relevance of neutrons 
to the 20 day bumps observed in a few GRBs, as the time coincides with i?traii/c 

Neutron signatures should be absent if the fireball is strongly dominated by a Poynting flux 
and has extremely low baryon loading. Then the neutron component decouples early, with a 
modest Lorentz factor r„, and decays at small radii. The upper bound on T n due to decoupling is 
r„ 300(Mq/10 26 ) 1 / 3 where Mn [g/s] is the mass outflow rate per unit solid angle of the fireball 
(see § 4.1 in Paper 1). 

We focused here on the neutron front and did not account for the 7-ray precursor that impacts 
the blast wave dynamics at R < i? acc = 0.7 x 10 16 (^ 7 /10 53 ) 1 / 2 cm, where £ 7 [erg] is the isotropic 
energy of the GRB (see Beloborodov 2002 and refs. therein). The analysis in this Letter is strictly 
valid for afterglows emitted at R > R& cc . Then the radiation- front effects, including the gap opening 
at R < 0.3i? acc , occur at smaller radii, and apply to an earlier afterglow. For a dense medium, 
where -Rdec < -Race, effects of the neutron and 7-ray fronts should be studied together. 

I am grateful to P. Goldreich, R. Blandford, and TAPIR group for their hospitality during my 
stay at Caltech, where part of this work was done. This research was supported by NSERC and 
RFBR grant 00-02-16135. 
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